Probing the ground state in gauge theories 
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We consider two very different models of the flux tube linking two heavy quarks: a string linking 
the matter fields and a Coulombic description of two separately gauge invariant charges. We compare 
how close they are to the unknown true ground state in compact U(l) and the SU(2) Higgs model. 
Simulations in compact U(l) show that the string description is better in the confined phase but 
the Coulombic description is best in the deconfined phase; the last result is shown to agree with 
analytical calculations. Surprisingly in the non-abelian theory the Coulombic description is better 
in both the Higgs and confined phases. This indicates a significant difference in the width of the 
flux tubes in the two theories. 



I. INTRODUCTION 

We can naively think of a heavy meson as being made 
up of a static quark and a static anti-quark separated 
by some distance. Around these fermions there will be 
a cloud of glue whose form is unknown but is typically 
thought of as being 'cigar shaped' (we neglect light quark 
flavours in this paper). 

It is important to realise that the matter fields in an 
interacting gauge theory which are not locally gauge in- 
variant, cannot be directly interpreted as observables. 
Thus a crucial role of the glue around quarks is to make 
the system gauge invariant. It is easy to imagine doing 
this by linking the fermions by a string like gluonic line. 
In an operator approach this is the familiar path ordered 
exponential which, assuming it to have a non-zero over- 
lap with the ground state, is evolved in time to form 
the rectangular Wilson loop approach to the interquark 
potential. 

It is worth noting that in such a string like state the 
only gauge invariant object is the overall colourless me- 
son state. There are no separately gauge invariant con- 
stituent quarks and, even at small interquark separations, 
there is no sign of a Coulombic potential. 

While producing a confining potential might be 
thought to compensate for this deficit, it is well known 
that lattice simulations are greatly improved by 'smear- 
ing' the string. This improvement is due to the naive 
Wilson loop being formed from an initial state where the 
flux is trapped on an infinitcsimally thin path. As we 
will see below in U(l) theory, where we can calculate an- 
alytically this leads to the stringy description actually 
having zero overlap with the true ground state unless a 
UV cutoff (such as a finite lattice spacing) is imposed. 
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An alternative description, which may be expected to 
be closer to the ground state at shorter separations in 
QCD, is to 'dress' the quarks separately so that they 
are each locally gauge invariant. For static quarks this 
may be done via a Coulombic dressing which has been 
shown [l[ to generate the perturbative interquark poten- 
tial (parts of which have been verified to next to next 
to leading order). It also allows an identification of 
the gluonic structures which generate the anti-screening 
and screening interactions, and incorporating this dress- 
ing removes perturbative infra-red divergences in on-shell 
Green's functions. Beyond perturbation theory the Gri- 
bov ambiguity sets a fundamental limit 2] on the observ- 
ability of individual dressed quarks, but it has previously 
[3] been shown that such Coulombic states may never- 
theless be used on the lattice simply by averaging over 
Gribov copies. 

In this paper we will compare the string like (or 'axial') 
and Coulombic descriptions of a heavy quark-antiquark 
system. We will do this analytically for deconfined U(l) 
and on a lattice for both compact U(l) and the SU(2) 
Higgs model. We will see that there are great differences 
in the two theories. 

We make contact between our analytic and numerical 
studies of the ground state through a particular probe, 
namely the ratio 



KV>|0}| 2 



(i) 



for any two states | yb ) and | ip' ) and where | ) is the 
(typically unknown) ground state of the theory. Simply, 
if this ratio is less than one then the overlap of the ground 
state with | yb ) is better than the overlap with | yb' ) , and 
we should expect | -0 ) to yield the closer description of 
physics in the ground state. If the ratio is greater than 
one then | yb' ) has the better overlap. 

This ratio can be realised both analytically for the 
deconfined U(l) theory and on the lattice through cal- 
culation of the large Euclidean time limit 



lim 

t — >oo 



(y/|e- gt |y/) 

(V>|e- At | yb) 



(2) 
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since in Euclidean space the large time limit is a vacuum 
projector 

This paper is organised as follows. We begin in Sec- 
tion |TT] by analysing a system of two non-confined abelian 
charges. This explicit investigation will allow us to build 
up ideas and methods which will be applicable more gen- 
erally. We will see that the Coulombic description yields 
the ground state while the axial system does not yield 
the correct potential and is severely divergent. This is in 
good agreement with physical expectations. In Section 
IIIII we describe the lattice techniques we employ and go 
on to examine the ground state in compact U(l) in both 
the deconfincd and confining phases. It will be shown 
that the simulations in the deconfined phase are in good 
agreement with the analytical results. In Section Hvl we 
apply our lattice techniques to an SU(2) Higgs theory. 
Finally in Section [V] we give our conclusions. 



gauge invariance by imposing Gauss' law on the Hilbert 
space, 

{AfaEjW %t) = _i0, *_*[.4,t] - p(z)9[A,t] , 

. (6) 

where p(z) is the charge distribution of the fermionic 
fields. 

We will frequently use the transverse (T) and longitu- 
dinal (L) projectors, which in momentum space are 

and obey T + L = 1, TL = 0. Transverse fields and 
derivatives are defined in a natural way by 

ffto-TaWfjfa), (8) 

and similarly for longitudinal parts. 



II. ANALYTIC TREATMENT OF 
DECONFINED U(l) THEORY 

In this section we consider QED with very heavy 
fcrmions, identified with static external sources. This 
is an ideal testing ground for our methods, as we may 
solve this theory exactly and analytically. This will al- 
low us to build up an intuition for the physics involved 
before proceeding to more complex systems. The iden- 
tification of the heavy fcrmions with physical objects is, 
however, somewhat suspect, as the Lagrangian fermions 
are not gauge invariant. Creating physical charges in 
gauge theories is the focus of the dressing approach, to 
be discussed below. We work in Weyl gauge, A° = 0, [|[ 
and with very heavy fermions our Hamiltonian is simply 
that of free gauge fields, 



H=±J d 3 x(E 2 + B 2 ) 



(3) 



where B is the magnetic field. Quantisation in Weyl 
gauge is straightforward 0, 0] as gauge potentials and 
electric fields are canonical variables. Hence we impose 
the non-zero equal time commutator 



[Ai{x,t),Ej{y,t)] 



iSij5 3 (x 



V) 



(4) 



We will work in the Schrodinger representation where 
the commutator ((U is realised on the state space by di- 
agonalising the field operator A^ and taking E^ to act 
as a derivative, in direct analogy to the position repre- 
sentation of quantum mechanics. States l^t) are then 
identified with functionals ^[A, t] of a vector field A on 
which the operators act as 



lA\A i (x)\*,t)=A i (xWAi\ 
S 



(A\Ei(x)\*,t) = -i 



SAi(x) 



*[At] 



(5) 



As we quantise the full gauge potential including its lon- 
gitudinal gauge non- invariant part, we must implement 



A. Dressed states in the abelian theory 

The Lagrangian fermions are not gauge invariant. In 
order to discuss physical, gauge invariant objects we dress 
the fermions with a function of the gauge field, which we 
write exp(Vt / [«4.]). We will be interested in states with 
two heavy fermions, where the state takes the form Q 



q(x 2 )q\x 1 )\0) 



(9) 



for q and q^ the heavy fermions and | ) the vacuum. In 
this section we will construct the most general functional 
W[.A] which makes the state invariant under the residual 
gauge transformations of Weyl gauge [47| , 



Aj(x) -> Aj(x) = Aj{x) + djX(x) 



q(x) 



e-' eX( - x) q{x), 



q\x) 



(10) 



(11) 



These imply that VF[«4.] must transform as 

W[A X ] = W[A] + ie \{x 2 ) - ie A(aci) . 
Equivalently, Gauss' law implies 

"JTFm W W = ieS 3 (x-x 2 )-ieS 3 (x- Xl ) , (12) 
ox 1 oAi(x) 

which constrains only that piece of W[«4] which is linear 
in A to have longitudinal part 



ie ^djAj(x 2 ) - ie —.djA^xx) 



(13) 



Here, Aj may be replaced with its longitudinal part 
Af. This is called the 'Coulomb dressing', first pro- 
posed by Dirac for modelling the gauge invariant static 
electron [ll|. The general dressing which gives a gauge 
invariant charge-anticharge state is therefore 



exp 



1 



ie ■ 



d j Aj(x 2 ) 



ie^djAjixr) 



J[A 7 



(14) 



for arbitrary J. We will investigate some of the choices 
of J[A T ] below. 
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B. The axial state 

What other principles shall we use to construct our 
dressings? Let us make an educated guess as to what a 
gauge invariant state might look like. Taking a bottom- 
up approach, we could think of beginning with the 
ground state in the vacuum sector, i.e. in the presence of 
no matter, p{z) = 0. This is the lowest, gauge invariant, 
eigenstate of our Hamiltonian ^ which is easily found 
to be given by the Gaussian 



of the string in the two directions perpendicular to n. 
The potential diverges as we remove the cutoff, and the 
axially dressed state is therefore infinitely excited [12]. 

In the classical theory it is known that the string be- 
tween the sources radiates energy and decays in time to 
an energetically more favourable state describing sources 
surrounded by a Coulomb field [l3[ ■ We will consider the 
quantum time development of this state below, but in 
preparation we first construct the true ground state. 
C. The ground state 



*o[^ T ] :=Det 1 / 4 v / ^V2 

i r d 3 P 



x exp 



(2tt) 3 



M-p)\p\T l3 (p)A J (p) 



(15) 



Gauss' law implying that this state is independent of Aj 
[48| . Now we add our fermion-antifermion pair, and ask 
what else we must add to make them gauge invariant. A 
common approach is to connect the sources by a string. 
This gives us what we call the 'axially dressed' state x> 



\[A] = exp 



ie I &z % Ai{z) 



(16) 



c 



We will take the curve C to be the straight line from x\ 
to x 2 . Fourier transforming, it is easily checked that x 
is of the general form (fT4|) , with J being the line integral 
of AF . Therefore x\>A\ q{x 2 )q i (x±)\ } describes a gauge 
invariant state. 

Have we found the ground state in the fermion- 
antifermion sector? Far from it - the axially dressed 
state (fT6|) is not even an eigenstate of the Hamiltonian, 
which, imposing a large momentum cutoff |p| < A, has 
expectation value 



sci 



(17) 



where the ellipses denote terms vanishing as A 
The first terms denote the vacuum energy, 



d 3 a 



d 3 p 



\P\ 



(18) 



which obviously represents the total energy of a noninter- 
acting photon 'gas' filling all space, and self-energies of 
the sources. The final term in (|17[) gives a linearly rising 
potential between the sources - i.e. they appear to be 
confined! The cutoff A 2 regulates a short distance diver- 
gence S 2 (0) coming from the infinitely small extension 



Gauss' law has already fixed the longitudinal piece of 
our wavefunctionals to be the Coulomb dressing, which 
is an eigenstate of the longitudinal part of the Hamil- 
tonian. The lowest eigenvalue of the transverse part of 
the Hamiltonian is just the vacuum functional \&o[>^-]- 
The ground state, or Coulomb dressed state, <E> in the 
fermion-antifermion sector is therefore 0, 0, EH 



<5>[A\ = exp 



ie 



1 



1 



djAj(x 2 )-ie ^d 3 Aj(xi) 



(19) 

with $[*4] q{x 2 )q\xi) describing a static, gauge invari- 
ant charge and anti-charge of minimal energy. The total 
energy of the ground state is found to be IR finite but 
UV divergent. For large values of the momentum cutoff 
we find 



E, 



c 



e 2 A 
2^ 



4ir\x 2 - cci 



i + e>(A- r 



(20) 



Again we have the vacuum and self energies. The fi- 
nal term is finite as the cutoff is removed and gives the 
Coulomb potential between the sources. The ground 
state in the fermion-antifermion sector is therefore de- 
scribed by two individually gauge invariant sources sur- 
rounded by a Coulomb field. 

Our dressed states ([TBI and (fT9"|) are related by 



X[A] = exp 



ie /dzMf(z) 



(21) 



It is this extra transverse piece which, in the classical 
theory, radiates away. We now study the time develop- 
ment of our quantum states. In our functional picture, 
time evolution described by the Schrodinger equation is 
implemented by taking the inner product of states with 
the Schrodinger functional. For the explicit form of this 
functional see, e.g. OElj]. Here, the calculation reduces 
to performing a Gaussian integral, with the result that 
the axial state at time t is 
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x[A,t] 
Af(t) 



-iE c t 



$[A] exp 



np 



nkAlip) 



d»p l re _ 2l|p|t ( l-cosp-fo-sO , niT .. (p)n . > 



(2tt)3 |p| 



(n-p) 2 



(22) 
(23) 



r 



which defines Af(t). There is no divergence at np = in 
any of our expressions, as the difference of exponentials, 
above, always ap pea rs to 'regulate' the notorious axial 
gauge divergence [17] . 

The axially dressed state, at time t, is comprised of 
the Coulomb dressed state together with an additional 
transverse contribution and a time dependent normalisa- 
tion j\f(t). Imposing a momentum cutoff, the transverse 
term and Af(t) are regulated expressions which are sup- 
pressed by rapid oscillations as t — > oo. Alternatively, we 
may use Abel's limit [II []J], 



OO 

lim f(t) = lim e / ds f(s) 

t—>oo e— 10+ J 



(24) 



of the Coulombic energy, we find (with a chosen regular- 
isation understood) 



lim c 

t — >oo 



iE c t 



x[A, t] 



_ d -JV(0)/2 



*[A] 



(26) 



This result echoes that of the classical theory. We have 
that both classically and quantum mechanically the ax- 
ially dressed state decays to the Coulomb state in the 
large time limit. The additional factor exp(— A/"(0)/2) 
ensures conservation of probability. 



D. Confined and deconfined U(l) charges 



to study the large time behaviour of these oscillatory ex- 
pressions. Physically, this limiting procedure adds a tiny 
imaginary part to the energies in the Fourier represen- 
tation of / which leads to exponential damping and the 
appearance of imaginary energy denominators. Apply- 
ing this, for example, to the transverse term we take the 
limit 



lim 



d 3 p —it 



(2tt)3 |p| 



ic 



np 



nkAlip) . 



This is very similar to our original axial dressing but with 
improved UV convergence. As e — > this expression van- 
ishes, and the transverse contribution, along with A/"(t), 
drops out. Therefore, pre-multiplying by the exponential 



We have seen that the lowest energy state of a gauge in- 
variant fcrmion-antifcrmion pair in U(l) gauge theory is 
given by the Coulomb dressing, the inter-fermion poten- 
tial being the expected Coulomb potential. Connecting 
the fermions by a string gives an unphysical, infinitely ex- 
cited state. It describes a confining potential between the 
matter sources and decays to the energetically favourable 
ground state. 

Having found the ground state exactly in this sim- 
ple theory, we may ask what use, if any, has the in- 
finitely excited axial state? To answer this, note that 
the Coulomb state $>[A]q(x2)q' (xi)\ 0) describes two in- 
dividually gauge invariant fermions. This may be seen 
from (11?)) , as the Coulomb dressing factorises into two 
parts, 



<5>[A] = exp 



exp 



^ [A J 



(27) 



The first exponential dresses q{x2), giving us a gauge 
invariant fermion, and the second exponential makes 
q*(xi) gauge invariant. In a deconfined theory this is 
the situation we expect — individual charges exist and 
are separately gauge invariant. 

It is not possible to unambiguously factorise the axial 
state into, similarly, one dressing for the charge and one 
for the anticharge. Taking l|16p. we Fourier transform 
and separate A into transverse and longitudinal parts, 



writing the line integral as 



d 3 p 
(2tt) 3 



inp 



VkA^ip) 



(28) 



n k Al{p) 



We have written out the longitudinal projector explic- 
itly. The first term is precisely the Coulomb dressing, 
as is necessary for gauge invariance. However, it is not 
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possible to factorise the second, transverse, part of (1281 
in a fashion similar to (|27p without introducing new di 
vergences. For example, we could try to write, 



inp 



inp inp 



(29) 



but whereas the left hand side of this expression is regular 
everywhere the right hand side requires the regularisation 
of axial gauge poles. 

What is the physical significance of this result? We 
saw earlier that the axial dressing appeared to lead to a 
confining potential. Physically, we therefore interpret the 
lack of a natural factorisation as an absence of individ- 
ually gauge invariant electrons and positrons. Instead, 
there is only a single, neutrally charged, object bound by 
a (thin) flux tube. It would seem, then, that in a confin- 
ing theory our axial state could be a better model of the 
true ground state than the Coulomb state. 

A suitable toy model of such a confining theory is 
compact QED, originally introduced by Polyakov [20( 
who argued that the Ad lattice model has a confinement- 
deconfinement phase transition. This was later proven 
analytically [M, HI HI, see also 0] and [H|. In the 
confining phase we expect the electric field to be concen- 
trated along a flux tube connecting the sources. Quali- 
tatively, this does not look much like the Coulombic field 
which we have seen is optimal so far, so the ground state 
could well differ significantly from the simple functional 

Although in this phase neither of our states will rep- 
resent the true ground state we may use them to extract 
information about it. We describe this method below 
and proceed to apply it to compact QED in the follow- 



ing section. 



E. A ground state probe 

What object are we to calculate which probes the 
ground state? Rotating to Euclidean space we consider, 
for any two states | ip ) and } , the ratio 



< Vle-^'l^) 



(30) 



In the large time limit contributions from excited states 
are exponentially suppressed. This limit then projects 
onto the ground state, call it | ), so that our ratio tends 
to 



l< -0 I 0>P 



(31) 



If the ratio is less than one in the limit then tp has the 
better overlap with the ground state, whereas if the ratio 
is greater than one ijj' is closer to the ground state. 

Below we will examine the ratio for ip' = x an d ip — 
our axial and Coulomb states respectively. Writing r = 
\x2 — %i\ the ratio of interest to us is 



R(r,t) 



{X\e- Ht \x) 
($|e-#*|$) 



(32) 



where the r dependence is contained in the states. In the 
deconfined U(l) theory we have studied so far is 



R(r,t) 



e Ect (x\e~ m \ X ) 

A 



exp 



d 3 p 1 



(2^)3 |p| 



(1 



,., , i 1 — cos(r np) . 



(33) 



exp(-AA(0)) as t 



oo , 



where the first line follows from <I> being an eigenstate. 

In the large time limit this correctly reproduces the 
mod squared overlap between the axial state and the 
ground (Coulomb) state, 

(x|$)=exp(-AA(0)/2) . (34) 

Af(0) is UV divergent (so that the overlap formally van- 
ishes, as one may expect for the overlap between a phys- 
ical and an unphysical state), but with the UV regulator 
in place Af(0) is positive definite and so our ratio be- 
comes | ( x I 1 2 < 1) as we would expect since | $ ) is the 
ground state and | x } is n °t- 



The ratio is plotted in Figure [T] To make contact 
with the lattice data of the following section, we take 
er 2 = 1.05 and measure distance in units of ir/A. In 
these dimensionless units momentum is cut off at \p\ < tt 
so that 1/A is associated with the lattice spacing. In- 
creasing the value of the cutoff does not qualitatively 
change the features of the plots. The integrand above 
is positive semi definite so that (with the regulator in 
place), we have R(r,t) < 1 for all r and t, as it must be 
since $ is the ground state and x is not. This may be seen 
in Figure [U which also shows that the ratio decreases 
rapidly with t to its asymptotic value of exp(— Af(0)). 
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FIG. 1: The ratio R(r,t) defined in ([32j) for unconfined 
U(l). Units are discussed below (|33p . As the separation of 
the charges increases the axial state becomes an increasingly 
poorer description of the ground state. 



We also observe that the ratio decreases as we increase 
the separation r of the charges. Increasing the length 
of the string concentrates more and more energy along 
the lengthening flux tube, giving an increasingly worse 
description of the true 1/r falloff of the Coulomb field. 
We turn now to the confining phase. 



III. A SIMPLE CONFINING MODEL: 
COMPACT U(l) 



where (3 = 1/e 2 with e the bare gauge coupling, 
plaquctte angles are defined by 



The 



0p(x)^ = + 9 v (x + n)- O^x + v)- 9 v (x) . (36) 

Because of the compact domain of support for the degrees 
of freedom, the U(l) theory admits magnetic monopoles. 
High precision measurements indicate that at 8 = /3 C rit, 
with /3 cr j t = 1.0111331(15) in 3 + 1 dimensions [2g|, there 
is a phase transition. Below the transition, Q < /3 C rit> 
there is an abundance of magnetic monopoles leading 
to confinement of electric charges by means of the dual 
Meissner effect [2(| ■ This has been convincingly con- 
firmed by lattice simulations [1^, |3(|. For 8 > /3 C rit, 
monopole nucleation is suppressed, and the theory is re- 
alised in the standard Coulomb phase. 

Our lattice comprises 12 4 points at /3=1.0 (confining 
phase) and /3=1.05 (Coulomb phase). A subtlety is that 
we must use open boundary conditions in the spatial di- 
rections and periodic boundary conditions in the tem- 
poral direction. This particular setup is necessary for a 
proper implementation of the axial gauge, further dis- 
cussed below. In order to reduce the influence of bound- 
ary effects, we always allow a distance of two lattice spac- 
ings to the boundaries when measuring observables. 

A standard heat-bath algorithm combined with micro- 
canonical reflections was used to bring the initial configu- 
rations into thermal equilibrium. Both disordered (hot) 
and ordered (cold) initial configurations were used and 
it was verified that our measurements were independent 
of the choice of the initial configuration. Measurements 
were taken with 20,000 configurations for each of the two 
phases. 



There is no exact analytic treatment of compact QED, 
although we note that the functional methods we em ploy 
have been successfully applied in 2+1 dimensions [2tll27l|. 
We therefore turn to lattice methods. In the following we 
outline our numerical approach and go on to calculate the 
ratio R{r, t). In the deconfincd phase we will sec that our 
analytic results can be reproduced numerically, providing 
a check on our methods. We will then go on to study the 
confining phase. 



A. Setting up compact QED on the lattice 

In the lattice formulation, the degrees of freedom are 
the fields 

Ufj,(x) = exp(i0 M (a:)) , —n < 6^(x) < tt , 

which are associated with the links of the 4-dimensional 
Euclidean space-time lattice. Using the standard Wilson 
action, the partition function is given by 

Z = J Ve» expj/^^cos^p^W)} » (35) 



B. Gauge fixing: Coulomb and axial gauges 

The quantities of interest are correlators of short 
Polyakov lines of temporal extent t and spatial separa- 
tion r, which begin in the time slice x° — and end in 
the time slice x° = t. Bringing these correlators into spe- 
cific gauges amounts to choosing specific dressings for the 
static test charge (and corresponding anti-charge). For 
more details see [3j. 

Let us focus first on Coulomb gauge fixing. Given the 
gauge transformation fl(x) — exp(ia(x)), 

Ufi{x) — > U^(x) = exp ( — ia(x + fi) + ia(x)) U^(x) , 
— 7r < a(x) < IT , 

(37) 

Coulomb gauge fixing is implemented by maximising, 
with respect to a, the gauge fixing functional 



T\U n ] = - Re 
1 J 2 



(38) 



X 2—1 

We use a standard iteration-overrelaxation scheme for 
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this task. To provide a stopping criterion for the itera- 
tion, we define 



A = — V A 2 



(39) 



A. r = Im^ r ^ 



U\ l {x) + U?\x-i)\ , 



i=l 



where the sum over (interior) lattice points excludes 
points on the boundary, and where consequently -/Vj n — 
(JVj — 2) 3 x N t . Introducing the photon field A^(x) via 



U M (x) = exp (ia Apix)) , 



(40) 



the above quantity A is a measure of the violation of the 
transversality condition, as can be seen from 



(41) 



We use A < 10~ 12 in our simulations. 

While Coulomb gauge fixing involves a non-linear opti- 
misation problem, axial gauge fixing can be implemented 
straightforwardly. The axial gauge condition is given by 



A 3 (x) = , equivalently U 3 (x) = 1 , 



(42) 



choosing our earlier unit vector n to point in the 3— 
direction. Let us assume that we have N 3 lattice points in 
the 3-direction, and the gauge transformations are num- 
bered according to 



flk = fi(xfc) with Xk = (x 1 , x 2 , ka, x 4 ) 



(43) 



The crucial observation is that, because of open bound- 
ary conditions, the gauge transformations fii and Qn 3 
are independent. Hence, it is always possible to choose 
£lk+i to satisfy 



U' 3 l (x k )=Sl k U 3 (x k )Sll +1 = 1. 



(44) 



As a result fl± remains undetermined and represents a 
residual gauge degree of freedom. 



C. Probing the ground state of compact QED 

Now that we have our lattice implementation of the 
dressings we may investigate their physical significance 
by looking at the overlap of the dressed states with the 
ground state. We will calculate the ratio R(r,t), defined 
in (|32ll . which probes the ground state in the large t limit. 
It is useful to introduce the matrix of transition ampli- 
tudes 



M{r, t) 



_ (C**(r,t) C x $(r,t) 



(45) 



where Cm^ :— {ijS \e \ip). This enables us to de- 
termine (for fixed values of r) the minimal value of t for 



which the contributions of excited states to the C$/$ (r, t) 
are sufficiently suppressed so that we can extract the 
overlaps with the ground state. Inserting a complete set 
of states, we obtain 



C^ V) =^e- £; »*(V''|n)(n|V) 



(46) 



where as usual E n denotes the eigenvalue of H for the 
state I n ). All contributions but that of the ground state 
(n = 0) to C^ij, can be neglected due to exponential 
suppression in the large t limit. A calculation of the 
determinant of the matrix M yields the large time limit 



a 2E Q t 



det[M(r,t)]- ,iWl l \ KX ,j U ^ 1 =0- 
/$|0)(0|x) l(x|o)| 2 

(47) 

Therefore, if the value of the determinant of M substan- 
tially deviates from zero (within error bars), the large t 
limit has not yet been reached. 

The calculation of C$5>(r,t), the Coulomb dressed 
Polyakov line correlator (pictorially, we may represent 
this object by |, see Q), is straightforward. Calculat- 
ing its axially dressed counterpart amounts to calculating 
the ordinary, gauge invariant, Wilson loop. This can be 
seen as follows: our starting point is the correlator of 
two Polyakov lines P(xi,X2) of temporal extent t and 
distance r calculated from gauge fixed links 



;x|o)(o|$> 

|0)| 2 



P(x u x 2 )[U 



X 2 = Xi 



rn 



(48) 



where n is the unit vector in 3-direction. Since in axial 
gauge the links U 3 are gauged to the unit elements, the 
above correlator can be completed to an r x t Wilson 
loop lying in the 3— t plane. Since the Wilson loop is 
manifestly gauge invariant, we find: 



p(x 1 ,xM] = w[u9] = w[u li 



(49) 



In other words, the loose ends of the Polyakov lines are 
joined by a straight line of parallel transporters U l _ l (x) in 
the initial and final time slices (this may be represented 
by □). Calculating the correlator of charges that evolve 
from axial gauge into Coulomb gauge, C$ x (r, t), means 
putting in a straight line of parallel transporters between 
the static charges in the initial time slice, the resulting 
object being of a U shape. Calculating the correlator 
from Coulomb to axially dressed charges C x $(r, t) pro- 
ceeds similarly, looking like a staple, IT 



D. Numerical results 

Calculating the familiar static potential between both 
axially and Coulomb dressed charges below and above 
/3 C rit, the phase transition is clearly illustrated. These 
results are displayed in Figure [5] (upper panel). We see 
the Coulomb potential above the transition, (3 = 1.05, 
and the linear rising potential below the phase transition, 
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FIG. 2: Upper panel: the U(l) inter-fermion potential in the 
confined and deconfined phases. Lower panel: behaviour of 
the determinant det[M(r, t)]. 



P = 1.0 (see also [3l| for numerical investigations of the 
confining phase) . In the lower panel of Figure [2] we have 
plotted the behaviour of the determinant det[M(r, t)]. 
We see that the determinant vanishes rapidly, coming 
effectively to zero at around t — 3a both above and below 
the phase transition. Thus the 'large time limit', in which 
we can extract good overlaps with the ground state, is 
reached quite rapidly. 

The numerical results for the overlap ratio R(r, t) are 
shown in Figure [3l The upper panel shows the Coulomb 
phase results (/? =1.05). Figure Q] is overlaid for com- 
parison. We see a close agreement between the ana- 
lytic solution and the lattice data, despite the differ- 
ent short distance cutoffs in use (the lattice spacing a 
in each direction compared to the 0(3) symmetric mo- 
mentum cutoff \p\ < A). We observe the formation of 
a plateau at around t — 3a, where the ratio approaches 
its asymptotic value, in agreement with the behaviour 
of the determinant of M (r, t) , cf . Figure O Increasing 



the spatial separation r of the charges probes the phys- 
ically more interesting long-distance behaviour. In the 
Coulomb phase, as seen in the analytic calculations, the 
ratio R(r, t) decreases as the separation r of the charges 
is increased. The energy of the axial state increases lin- 
early with r, giving an increasingly worse overlap with 
the ground (Coulomb) state. 

Encouraged by the agreement between our numerical 
and analytic results, we now consider the confining phase, 
j3 = 1.0. The results are shown in the lower panel of Fig- 
ure [31 Most notably, the ratio is now greater than one, 
which signals that the axially dressed state has the bet- 
ter overlap with the ground state. Although the signal 
is rather noisy in this phase, compared to the Coulomb 
phase, one can still observe the formation of a plateau 
in R{r, t) up to r = 3a. Beyond this value, the statistics 
were not sufficient to extract data points which met our 
acceptance criterion of a relative error smaller than 10%. 
As r increases the ratio R rises well above 1. There is no 
string breaking, and for larger r the Coulomb potential, 
falling off as 1/r, becomes a poorer and poorer descrip- 
tion of the long flux tube connecting the sources. These 
results show that the ground states differ significantly 
above and below the phase transition. Our physical in- 
tuition tells us that in a deconfined theory we expect to 
find two separate, gauge invariant charges. In a confining 
theory, however, we presume that physical states should 
be composite, uncharged objects. These expectations are 
indeed supported by our results. In the confining phase 
the Coulomb state, with its individual physical charges, 
is not a good description of the ground state. Instead it 
is now the, overall uncharged, axial state with its nar- 
row flux tube (its divergences regulated by the lattice 
spacing) which more closely resembles the true physics. 

In the next section we will extend our discussion to 
non-abelian theories. 



IV. TOWARDS QCD: THE GROUND STATE IN 
THE SU(2) HIGGS MODEL 

Ideally, we would like to study the overlaps of the 
Coulomb and axially dressed states with the true ground 
state in a non-abelian pure gauge theory (no matter) 
which has both a confining and a non-confining phase. 
However, for an SU(N) pure gauge theory asymptotic 
freedom implies that there is no fixed point in coupling 
constant space other than that at (3 — * oo. Hence, the 
SU(N) pure gauge theory has only a single, confining, 
phase at zero temperature. High temperature SU(N) 
Yang-Mills theory should provide a non-confining frame- 
work. However, as the critical temperature of, e.g., SU(2) 
pure gauge theory is approximately T c « 300 MeV, the 
temporal extent of the lattice is 

L t w n ^] r , T « ~fm. (50) 



300 MeV 3 
Hence, excited states are only suppressed by a factor 
exp(— H L t ) which is not enough to provide a good over- 
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FIG. 3: The overlap ratio R(r,t) as a function of t and r 
in the deconfined/Coulomb phase, above the phase transition 
(upper panel, solid lines are analytic results) and in the confin- 
ing phase, below the phase transition (lower panel). Lattice: 
12 4 , spatially open boundary conditions, temporally periodic 
boundary conditions, /3 — 1.05 (upper panel) and /3 = f.O 
(lower panel). 



lap with only the ground state. 

The simplest theory which provides both a confining 
and a non-confining phase at zero temperature is an 
SU(2) gauge theory with a Higgs field <fi in the fundamen- 
tal representation. This will be our model, but before we 
can study it we must briefly discuss what we mean by 
the Coulomb and axial dressings in a non-abelian theory. 



However, in a perturbative expansion the lowest order 
term of the non-abelian dressings is typically given by 
taking the dressings of Section [ill and replacing A.j with 
A^T a for T a the Lie algebra generators. For example, 
the Coulomb dressing is, to lowest order in the coupling 



exp g 



djAjT a 



(51) 



of similar form to the abelian case, (|T4|) . All orders in 
this dressing are required to make gauge invariant states, 
and for the construction of higher order terms, and issues 
relating to Gribov copies, we refer the reader to [1, H, H3] . 

Just as in the abelian case, correlators between dressed 
states are calculated (using the full dressing) by putting 
time slices into particular gauges - Coulomb and axial in 
our case. The methods employed in the U(l) theory may 
then equally be applied to non-abelian theories. With 
this understanding, we proceed to define and examine 
our SU(2) model. 



B. The SU(2) Higgs model 

We will study a gauged 0(2) model where the length 
of the Higgs field is restricted to unity, 



= 1 



(52) 



It has long been known that the confinement phase and 
the Higgs phase are smoothly connected [33| . A true or- 
der parameter, i.e., a local operator which distinguishes 
between both phases, does not exist. When static quark 
and antiquark sources are exposed to the gluon field the 
electric flux tube between them may break if the energy 
stored in the tube exceeds twice the mass of the Higgs 
particle, and the string tension vanishes. Although a lo- 
cal order parameter does not exist, the physics in the 
string breaking phase is vastly different from that in the 
Higgs phase: while short electric flux tubes form in the 
confining phase (thereby giving rise to a "string tension" 
at intermediate distances), there is no flux tube at all 
in the Higgs phase. Although this behaviour cannot be 
detected with local order parameters, non-local quanti- 
ties, such as the vortex percolation probability [34j |. can 
map out the phase diagram. In fact the critical line for 
vortex percolation follows the first order line of the phase 
diagram, and turns into a Kertesz line in the crossover 
regime [MS 111- 
The Higgs doublet 



A. Non-abelian dressings 

As in Section [Til we are interested in the properties 
of quark-antiquark states made gauge invariant by us- 
ing dressings. The form of the non-abelian dressings 
is much more complex than their abelian counterparts. 



h,2 e 



may be cquivalently written as a matrix a, 



oi -<p 2 



(53) 



(54) 
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where a is unitary. The lattice action of the theory is 
then 

S = Swilson + ^Higgs , (55) 

Swiison = /3 ^ tr U M (x) U v (x + fi) U^x + v) U%(z) , 

(56) 

S , Higgs = K^ ^ tr a ( x ) U^(x) a" 1 (x + fi) , (57) 



.r . / / 



where K is the Higgs hopping parameter. The "gluonic" 
degrees of freedom are represented by the SU(2) link ma- 
trices Ufj,(x) as usual. The action and partition function 
Z, 



Z = J VU^ Va exp (S) , 



(58) 



are invariant under SU(2) gauge transformations: 

u"(x) = n(x) u^(x) + , lie su{2) (59) 



! (x) = a(x) Stf(x) 



(60) 



In addition the theory possesses a global 0(2) symmetry: 

a"(x) =wa(x) , we 0(2). (61) 

The model may be studied using standard lattice tech- 
niques. We work with (3 = 2.2, several values of k and 
values of the lattice size N [49j . The susceptibility of the 
Higgs action, 



CHi, 



1 

TV? 



(62) 



may be used to explore the phase structure (for fixed 
[3). Our numerical results are shown in the upper panel 
of Figure |4j The susceptibility is rather independent of 
the lattice size and strongly peaked at a critical value 
k = k c , where k c rs 0.839(2). This indicates that the 
model undergoes a crossover from the confinement phase 
to the Higgs phase at k c . 

Figure 01 lower panel, shows the static quark poten- 
tial for k = 0.825 (confined phase) and for k = 0.88 
(Higgs phase) which is significantly bigger than the crit- 
ical value k c . In the low k phase, we observe a linear 
rise of the potential at large quark-antiquark distances. 
String breaking, induced by the presence of the funda- 
mental Higgs field, is not observed for ra < 8. We do, 
however, observe a string tension 

a a 2 w 0.15(1) , ,9 = 2.2, k = 0.825 , (63) 

that is somewhat reduced compared to the string tension 
of pure SU(2) Yang-Mills theory, aa 2 w 0.28(1) (for (3 = 
2.2, k = 0). In contrast, the potential data at k = 0.88 
are well fitted by a Coulomb law. 
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FIG. 4: Upper panel: The susceptibility of the Higgs action 
(|62p as a function of k for /3 — 2.2 for several lattice sizes N. 
Lower panel: The static quark potential for a 16 4 lattice and 
13 = 2.2 for two different k values. 



C. Overlaps with the ground state 

The ratio R{r, t) in this theory is shown in Figure[SJ In 
the Higgs phase (k — 0.88) the ratio is below one show- 
ing that the Coulomb state has a better overlap than the 
axial state, similar to our U(l) results. Comparing the 
two panels of Figure \E\ one sees that, for a given separa- 
tion r, the ratio is slightly larger in the confined phase 
(k = 0.825) but it is still below one. This implies that, 
counter intuitively, the Coulomb dressed quarks still pro- 
vide a better approximation to the true ground state. We 
attribute this finding to a bad overlap between the axial 
and ground states caused by the infinitely thin flux tube 
of the axial state. The true ground state is expected to 
sustain a flux tube with a thickness of the order of several 
lattice spacings. In this regime, the ansatz of a very thin 
flux tube, whose potential energy is only finite because 
of the lattice spacing, is not as good a description as the 
Coulomb dressing. 
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FIG. 5: Ratio of overlaps between the axially dressed and the 
Coulomb dressed quark states in the Higgs phase (k = 0.88, 
upper panel) and the confining phase (k = 0.825, lower panel). 
Lines have been added between data points to guide the eye. 



Our analysis suggests, then, that the traditional view 
of a thick flux tube [38| is a better description of the true 
ground state. In Figure [5] we see that as the separation 
of the charges increases, the ratio R decreases (at all 
times plotted), indicating that a thin string is a poorer 
description of the ground state at large r. To emphasise 
this point we plot, in Figure [51 our ratio in the 'large 
time limit', t = Aa, where the ratio is clearly seen to 
decrease with increasing separation in both phases. This 
indicates, in particular, that in the confining phase, at 
fixed lattice spacing, the thin string is a poor description 
of the ground state at large separations. 



V. CONCLUSIONS 

In this paper we have studied various locally gauge 
invariant ansatze for the lowest energy state in the 




FIG. 6: Ratio of overlaps at t — 4a, plotted as a function of 
the quark separation r/a in the Higgs phase (k = 0.88) and 
the confining phase (« = 0.825). 



heavy fermion-antifermion sector of both abelian and 
non-abelian gauge theories. One of these constructions 
corresponds to a Coulombic dressing around each matter 
field, i.e., the state is made up of two separate physi- 
cal charges (with opposite signs). In the other, axial, 
ansatz the fermions are linked by a string like flux tube. 
This dressing does not factorise naturally, so no locally 
invariant charges can be defined. There is instead a sin- 
gle, uncharged object. We have compared how these two 
constructions overlap with the ground state for different 
theories and in various phases. Our primary tool here 
is the ratio R(r, t) , defined in (p{2")) , which compares the 
overlaps between two ansatze and the true ground state. 

Firstly we considered U(l) theory where we have a 
good understanding and expect the ground state to cor- 
respond to two Coulomb charges. Using the functional 
Schrodinger representation we were able to construct the 
Coulomb state and also to show that the axial state is 
unstable and decays into the Coulomb one. This was fol- 
lowed by lattice simulations of compact U(l). In the de- 
confined phase the results agreed very well with our ana- 
lytic arguments, the Coulomb dressed state being clearly 
preferred. 

In the confined phase of compact U(l), where we do 
not have a good analytic understanding, our simulations 
showed that the axial state was preferred. This is per- 
haps what one might expect, the axial state corresponds 
to a flux tube between the fermions. However, it is an 
extremely thin flux tube — which is only non-zero on the 
links between the fermions — and in previous work Q 
it was shown that for SU(2) a Coulombic dressed state 
more readily yielded the interquark potential than an 
axial state. This seems to indicate a major difference 
between the confining flux tubes in compact U(l) and 
SU(2). 
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Spurred on by this we then studied the SU(2) Higgs 
model where there is a crossover from a Higgs phase to 
a confined phase. We calculated the potential on either 
side of the crossover and clearly saw the two phases and 
the effect of the Higgs field on the potential compared 
to the pure SU(2) theory. As with pure SU(2), 3, our 
simulations clearly demonstrated that the Coulomb state 
was preferred to the axial ansatz on both sides of the 
crossover. We interpret this as a consequence of the ex- 
treme thinness of the axial flux tube. Presumably the 
flux tube in the true ground state in this non-abelian 
theory is much thicker than its compact U(l) counter- 
part and this is why the Coulomb dressing is preferred 
despite its slow (1/r) fall-off. It is interesting to speculate 
on the phenomcnological implications of this insight. 

Our results agree with those of the investigations 
[3^ , |40| , where thickening of the tube with an increase in 
charge separation was observed at a fixed lattice spacing. 
However, those authors, based on an ansatz for the form 
of the string, also claimed that for finer lattices the width 
of the flux tube was proportional to the lattice spacing. 
This surprising result was argued to imply that in the 
continuum limit the width of the flux tube vanishes, i.e., 
the continuum ground state would be described by an 
infinitely thin string of flux. This result is the subject 
of some current debate and we will return to it in future 
work. 

The differences between the overlaps in the confining 



phases of the abelian and non-abelian theories clearly de- 
serve further study. It would be interesting to try to con- 
struct different dressings which correspond to a thicker 
flux tube ('cigar shape') around the fermions but which 
still fall off more rapidly than the Coulombic dressing. 
One possible way to approach this would be to use a 
dressing which corresponds to some form of interpola ting 
gauge, somewhat analogous to ideas of 't Hooft [4l|, |42| 
and others [H, 0, Hf| [40| on gauges depending on some 
'flow' parameter. A gauge interpolating between axial 
and Coulomb gauge, say, is expected to yield, for an op- 
timal choice of the interpolating parameter, a still better 
overlap with the ground state. It could then be used for 
more efficient lattice calculations of, e.g., the interquark 
potential. 

It would also be interesting to compare the overlap of 
axial and Coulombic dressings in the SU(3) theory; this 
would let us investigate whether modifying the number 
of colours significantly alters the flux tube between two 
heavy quarks. 
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